39045961-644d-4473-8873-08c220fedbad 21 3 5 auto

&[DATE] &[TIME] - &[FILENAME]

x f x x β 3 el x * β 2 el + β 4 el x 2^* + (* β 1 el + (* 1 x x β 7 el x * β 6 el + β 8 el x 2^* + (* β 5 el + (* + / : x Deriv x f (x diff : 21 line

β 4.69292878185 2.485645088668 0.005049502793 - 0.008399478307 - 5.19337552684 0.698811411628 0.043875966805 - 0.000673223801 - 81 mat :

u Q x f : z Int u Q u 0 x int : x Int z Int : 31 line

recast integ'r

t0 1 time :

Smath plot span => 164 pixels Integrator runs 164 times

0 x ≤ ( x 15 ≤ ( & x Int not sanity fit if 15 64 15 32 2 2 mat 2 1 sys

x c 10 x ≤ (x 15 ≤ (& cases :

1 time t0 - 2 s *

t0 1 time :

x f x c * x Deriv x c * 2 1 sys

1 time t0 - 0.7 s *

x Int @ integrator acc 50 ... timing 2/164 = 0.012 sec ≡

f a b n Simpson h b a - n / : 13 / h 2 / a f b f + (* h a k h * + f k 1 n 1 - sum * + (* 23 / h * a h 2 / - k h * + f k 1 n sum * + eval 21 line :

==============

===================

These two components in // ~ ½ time

Explode the derivative components

x Denom 1 x x x β 8 el x * β 7 el + (* β 6 el + (* β 5 el + (* + (2^:

t0 1 time :

x Denom x c *

1 time t0 - 14 / s *

x Num x x β 4 el x * β 3 el + x β 4 el x diff x * 3 β 4 el * + β 3 el x diff + (* + 2 β 3 el * + β 2 el x diff + (* 2 β 2 el * + β 1 el x diff + (* β 1 el + (1 x x x β 8 el x * β 7 el + (* β 6 el + (* β 5 el + (* + (* x x x β 4 el x * β 3 el + (* β 2 el + (* β 1 el + (* x x β 8 el x * β 7 el + x β 8 el x diff x * 3 β 8 el * + β 7 el x diff + (* + 2 β 7 el * + β 6 el x diff + (* 2 β 6 el * + β 5 el x diff + (* β 5 el + (* - :

t0 1 time :

x Num x c *

1 time t0 - 0.55 s *

Num x Num x Denom / x c * 2 1 sys